Parity-dependent proximity effect 
in superconductor/antiferromagnet heterostructures 
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We report the effect on the superconducting transition temperature (T c ) of a Nb film proximity coupled to the 
synthetic antiferromagnet Fe/{Cr/Fe}jv_i. We find that there is a parity dependance of T c on the total number 
of Fe layers, N; locally T c is always a maximum when N is even, and a minimum when N is odd. The Fe 
electron mean free path and coherence length are indicative of dirty limit behavior; as such, we numerically 
model our data using the linearized Usadel equations with good correlation. 

PACS numbers: 74.45.+C, 74.62.-c, 74.78.Db 



During the past fifteen years the coupling of thin film su- 



iting 

perconductors (S) to ferromagnets (FM) JjJ] has been exten- 
sively studied. In the most simplistic sense, these electron 
coupling phenomena can be considered to be mutually exclu- 
sive with spin alignment enforced by ferromagnetic exchange 
and anti-parallel spin alignment necessary for singlet BCS su- 
perconductivity. According to theory, the penetration depth of 
the superconducting proximity effect (SPE) into normal met- 
als (NM) is governed by the electron phase breaking length. 
When NM is substituted with a FM, the singlet superconduct- 
ing correlations and phase coherent effects are destroyed by 
the exchange interaction / within the FM coherence length 
£/ oc y/l/I because of the splitting of spin-up and spin-down 
conduction bands. For S/FM bilayers, a consequence of I on 
S is the non-monotonic dependence of T c on FM layer thick- 
ness (df) for when the superconductor is thinner than its BCS 
coherence length [Hfl. This can be understood in terms of an 
interference effect between the transmitted singlet pair wave 
function through the S/FM interface with the wave reflected 
at the opposite FM interface. 

Recently, Andersen et al. @] predicted a SPE dependence 
on the number of antiferromagnetic (AF) atomic planes (N) 
in S/AF/S junctions such that the junction's ground state is 
or 7r depending on whether N is even or odd. Experimen- 
tally, this SPE dependence on N is difficult to realize because 
it demands atomic thickness control of the AF. So far, exper- 
iments looking at the proximity of S to AFs have focused on 
thicker AF layers |01 where this parity dependence does not 
exist. A way to investigate a parity dependent SPE is to use 
synthetic antiferromagnets (SAFs) which exploit the AF cou- 
pling of FM layers separated by a NM spacer [JsJ] . 

So far, the proximity effect of SAFs coupled to S mate- 
rials have not been considered although the T c dependence 
of complicated multilayers have been extensively studied |@]; 
for example, the control of T c in pseudo-spin valve FM/S/FM 
structures has been proposed Q] and realized experimentally 
with mK differences in T c between parallel (P) and anti- 
parallel (AP) FM configurations. The most similar structure 
to the one we study here was proposed by Oh et al. J^] in 
which the S layer is in proximity to decoupled FM/NM/FM. 
In this Letter, we report the proximity-suppression of T c of a 
Nb film coupled to a Fe/{Cr/Fe}jy-i SAF. In doing this, we 
find that T c has a pronounced parity dependence on N such 



that when N is even, T c is a local maximum and, conversely, 
when N is odd, T c is a local minimum. 

Films were grown in Ar (1.5 Pa) on oxidized (120 nm) 
Si (100) (surface area: 5-10x5 mm 2 ) in a diffusion-pumped 
ultrahigh vacuum sputter deposition system, consisting of 
the following: three dc magnetrons; a computer operated 
(rotating) sample table; a liquid N2 cooling jacket; and a 
residual gas analyzer. The vacuum system was baked-out 
overnight prior to each experiment, reaching a base pressure 
of l-4xl0~ 6 Pa. One hour before depositing, the system 
was cooled with liquid N2 giving a final base pressure of 1- 
3xl0 -8 Pa, a residual O2 pressure of <4xl0 -9 Pa, and an 
outgassing rate of « 1 x 10~ 8 Pa s _1 . Targets (Nb, Fe, and Cr 
of 99.9% purity), were pre-sputtered for ~ 15 minutes to re- 
move contaminants from their surfaces and to further reduce 
the base pressure of the vacuum chamber by getter sputter- 
ing. To control film thickness (d) and to ensure clean inter- 
faces, films were grown in a single sweep by rotating the sub- 
strates around the symmetry of the chamber under stationary 
targets. Growth rates were pre-calibrated by growing films on 
patterned substrates and, from a lift-off step-edge, d for each 
material was established using an atomic force microscope. 
Typical growth rates with a sweep speed of 1 rpm per pass 
were: 1.0 nm for Cr; 0.6 nm for Fe; and 1.4 nm for Nb. The 
average film roughness was < 3 A over 1 /zm. 
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FIG. 1: (color online). Dependence of H aat on dc r in 
Nb(20nm)/Fe/{Cr/Fe} X 7 films for: d Fe =0.6 nm (□); and d F e=l.O 
nm (■). Solid lines are to guide the eye. Inset: illustration of 
S/FM/{NM/FM}]v-i structure with N FM layers and N - 1 NM 
layers. Vertical arrows indicate FM polarization direction. 
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The antiferromagnetic interlayer exchange coupling (AFC) 
between the Fe layers was investigated in two sets of 
Nb/Fe/{Cr/Fe} X 7 films: for dp e =0.6 nm; and for dp e =l.O 
nm. (ijvb was constant (20±1 nm) while dc r varied in the 
0.5-2.5 nm range. A 2.8 nm capping layer of Nb was grown 
on the top Fe layer. To quantify AFC in these films and to 
optimize dc r to give the largest AFC energy J (given by 
— 4J = noH sa tMdF e where M is the magnetization of the 
Fe), the required saturating field H sat needed to align the Fe 
layers was measured with a vibrating sample magnetometer at 
room temperature as a function of dc r ', see FigJTJ For dp e = 
0.6 nm, H sat is a maximum value of (1.16±0.04)xl0 6 A/m 
for d C r = 1.0±0.2 nm with M w 0.8xl0 6 A/m (similar to 
Hot]). This corresponds to - J w 1.8 xl0~ 4 Jm~ 2 , which 
compares well to previously reported energies Jit]. To achieve 
the most efficient AFC, we chose dp e = 0.6 nm and dcr = 
0.9 nm for the appropriate layer thicknesses. 

To establish a dMb range that is strongly affected by the 
presence of the SAF film, the T c of both bare Nb and 
Nb/Fe/{Cr/Fe} X 7 films was measured in the 15-45 nm range; 
see Fig. |2] The critical d^b is ~ 12 nm. T c was estimated by 
measuring the resistance of a film as a function of temperature 
(R(T)) using a standard four contact technique; for this, two 
instruments were used: a custom made liquid He dip probe 
and a pumped He-4 temperature insert. To make electrical 
contacts, films were ultrasonically wire-bonded onto copper 
carriers with Al wire (25 fim diameter). Samples were fixed 
to their carriers with silver conducting paint. An ac current of 
±10/iA was applied. T c was estimated from warming curve 
data and defined as the mid-point of the R(T) transition. 

Fig. Oshows the dependence of T c on N for two Nb thick- 
nesses: (a) 22 nm and (b) 18 nm. We see that T c is a local 
maximum when N is even and vice-versa. The magnitude of 
\T C (N + 1) — T C (N)\ decreases as N increases: the largest 
for T c (2) - T c (3) =1.38 K in (a) and 2.06 K in (b). Assum- 
ing dirty limit behaviour in the multilayers (£ > £), we have 
adap ted the Usadel equations, as discussed by Fominov et al. 
Il ill and Oh et al. JjJ, to model the situation; as illustrated in 
the inset of Fig. Q] 
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FIG. 2: (color online). Dependence of T c on djvt: bare Nb (A); two 
Nb/Fe/{Cr/Fe}x7 sets for dFe = 0.6 nm and dc T = 0.9 nm (□ and 
■); theory fit (•) explained in text. 
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FIG. 3: (color online). Experimental (■) and theoretical (□) depen- 
dence of T c on N: (a) djvi, = 22 nm; (b) djvfc — 18 nm. For com- 
parison, data points (0) are for a theoretical S/FM bilayer and (A) 
for bare Nb taken from Fig. 2. Inset: (i) the theoretical non-collinear 
case for various 29 angles between FM layers; (ii) illustration defin- 
ing 29 between two FM layers; (iii) typical resistivity curves for N = 
2,3 and 4 in the function of T. (c) Spread in T c (AT C ) defined as 
20%-80% of the transition. 



The normal Green function is G = sgn(wi) when T « T c 
and the linearized Usadel equations for the anomalous Green 
function F are: 

hD. d 2 F, 



2 dx 2 



h\uji\F s + A = 



(1) 



inside the S layer, 

HD f d 2 F fk 



h\wi\F fk - i(-l)*JBgn(w,)-P>ik = (2) 



2 dx 2 

inside the fcth (k = 1, 2, ...TV) FM layer and 

hD n d 2 F nk . . 

K\wi\F n k = 



(3) 



2 dx 2 

inside the fcth (k = 1, 2, ...N - 1) NM layer. 

The appropriate diffusivities are D s , Df and D n , while A 
is the superconducting pairing potential in the S layer. The 
(— l) fe factor before / accounts for the alternating polarization 
of the FM layers. The Usadel equations are supplemented by 
the self-consistency equation 
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H\uji 



- F, 



(4) 
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where the summation goes over the Matsubara frequencies 
lui = (21 + l)-Kk B T/h with all integers I. T is the T c of 
bare S of thickness d s , i.e. see Figj2] The anomalous function 
F obeys the boundary conditions 



dF a {Xa,b) dF b (x a ,b) F b (x a ,b) - F a {Xa,b) 



dx 



Ob- 



dx 



R a ,< 



(5) 



at the points x = x a< b separating any a and b layers; a and b 
can be either of the S, FM or NM layers and in each case we 
use the appropriate normal-state conductivities {a a and a b ) as 
well as the appropriate interfacial resistance per unit area R a b 
between a and b layers. Since all the FM and NM layers are 
identical, the three conductivities are <r s , <jf and er„ and the 
two interfacial resistances are R s j and i?/, n . 

If we introduce the formal vector $> b = (F b , dFb/dx), we 
can turn Eq. <J2J and ([3]i into vector equations for 3>b (b refers 
to an arbitrary layer again), which take the common form 



— * fc = M 5 * 5 

dx 

where the values of k? are 



with Mb = 



1 

k 2 



k 2 - 
K fk - 



2h\u l \+2i{~l) k Isgn{u Jl ) 

hD f 



and k 2 nk 



2N 

Dn 



(6) 



(7) 



in the fcth FM and the fcth NM layers, respectively. Eq. © 
is a simple linear differential equation. If is known at the 
left side (x = x a> b) of a layer, its value at the right side {x = 
x a ,b + db) can be expressed in terms of matrix exponentials 

&b(x a ,b+db) = A b &b{x a ,b) with A h = exp(M b d b ). (8) 

The boundary conditions (Eq. [5]) between any two layers a 
and b can be written in the vector form as 



®b(x a ,b) = B b , a $ a (x a , b ) with B b 



1 Ra,bOa 
(J a j Ob 



(9) 

The derivative of the anomalous function F vanishes at the left 
end of the multilayer, therefore *& fj^(—dfN — d n (N — 1)) = 
(C, 0), where C is a complex number. The formal vector 3? 
at the left side of the S layer is now expressed in the function 
of C by systematically going through the layers and using the 
appropriate matrices appearing in Eqs. (0 and (0: 



* a (0)=L(wO 



where 



(10) 



7V-1 



L(wj) = B fl) /i JJ [A/fcB/fc )n fcA n fcB nfe j( fe+1 )]A/jv (11) 

k=l 

is a complex matrix only dependent on material parameters 
and the Matsubara frequency oji, which depends on T. 

To obtain T c , we apply the multimode method developed 
by Fominov et al. Il ill . If we write down the components of 
the vector Eq. dTOb and eliminate C, we obtain 

dFM - Qi^FM with 0(^ = ^4, (12) 



which is analogous to Eq. (8) in 11 ill . From here, the prob- 
lem is reduced to finding the highest T for which the de- 
terminant of a matrix K is zero. If we include M modes 
(m = 1,2, ...Al) in addition to the single-mode method and 
take the first M + 1 Matsubara frequencies (I = 0,1, ...M), 
the elements of K are 



R(u>i) cos(fc d s ) - k a sin(k d s ) 

Am = : ^ and 

h(ui + k 2 D s /2) 



(13) 



_ R(uji) cosh(fc m rf s ) + fc m sinh(fc m 4) 

^-fc^ s /2) ' (14) 

where R(u>{) is given by 

[Im(QM)] 2 



R{loi) = Hfi(Q(w,)) 



Re(Q(wj)) + fco tanh(fe s d s ) 



(15) 



in the function of Q(uji) and fc s = ^j2ujijD s . The quantities 
fco and fc m (m = 1,2, ...M) are the smallest positive roots of 



T 



'I 5 



Ank B T 



(17) 



(Z.r 



which are obtained from the self-consistency Eq. and con- 
tain the digamma function ijj. 

Using this method, we determine T c numerically; the cal- 
culation is repeated with various values of T, for < T < Tq. 
By obtaining the roots fco an d k m of Eqs. ( fToT ) and ( fTTI i. and 
evaluating R(oji) through Eq. (15[ , the matrix K can be found 
for all T. The largest value of T for which dct(K) = corre- 
sponds to T c . The multimode method with Al — > oo is exact 
but in most cases the inclusion of M = 8 modes suffices. 

To apply this numerical method with the minimum num- 
ber of adjustable parameters, we have measured a for Nb, 
Fe and Cr thin films for T <10 K to be 1.9xl0 6 (fim)" 1 , 
6.6xl0 6 (fim) -1 , and2.2xl0 6 (f2m) -1 , respectively. From 
these values, we estimate the electron mean free paths of Fe 
and Cr via I = amvF-(ne 2 )~ l with vp the Fermi velocity; 
vp for Fe lfl2ll and Nb lfl3ll is taken from literature, while 
for Cr we assume a similar vp to Fe. The density number 
of electrons is n — (8ir/3).(mvp/h) 3 with m the electron 
mass, giving an I of 2.7 nm and 0.9 nm for Fe and Cr, re- 
spectively. For Nb, we determine I by choosing a value that 
gives the best fit (2.2 nm). With I known, we calculate D for 
Nb, Fe, and Cr via D = ivp/3, giving 2.2 x 10~ 4 m 2 s- 1 for 
Nb, 18xl0" 4 m 2 s- 1 for Fe, and 6.0xl0" 4 m 2 .^ 1 for Cr. 
Finally, we estimate the coherence lengths of Nb and Cr with 
£ s n = y/hD s , n /2irkBTo, giving 5.8 nm for Nb (similar to 
H) and 9.6 nm for Cr, while for Fe £/ = x /hDf/I giv- 
ing 3.7 nm (similar to ifioll ) assuming / ~ 1000 K fl5ll . For 
both Nb and Fe, £ > £, which is indicative of dirty limit be- 
havior and, therefore, justifies our use of the linearized Usadel 
equations. For the interfacial resistances we take Rp e ,Nb ~ 
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10~ 15 ilm 2 and Rf & ,Ct ~ 10~ 17 ttm 2 (the quality of the fit 
does not depend strongly on these parameters). The important 
material parameters used/calculated here are listed in Table J] 
With these values, the numerical model agrees well with the 
experimental data; see Fig. [3] 

We have generalized the model to consider the behavior 
when the Fe layers are non-collinear; we use the linearized 
Usadel equations containing both the singlet and triplet com- 
ponents of the anomalous Green function Il6ll . We find that 
with a small change in 29, the angle between the polarizations 
of adjacent Fe layers, the T c is more strongly reduced and the 
parity dependent oscillations fade away faster; see Fig. |3a,i). 
In the limiting case of 29 = n we recover the SAF behavior. 
An applied magnetic field will produce a non-collinear con- 
figuration; however, in our samples this could not be achieved 
without directly suppressing the Nb T c . With improved con- 
trol of J the field required to reorient the layers could be sub- 
stantially reduced. 

In conjunction with the experimental results, this demon- 
strates that a large change in T c can be obtained by switching 
the SAF between P and AP configurations. In the particular 
case of TV = 2, which corresponds to the Oh et al. Jgt] spin 
valve, the change in T c from AP to P is > IK, which is a 
many times higher than the changes experimentally observed 
in the analogous superconductor PS V structures 18fl. 

This Letter has shown that the parity of AFs with perfect 
order have a profound effect on the proximity effect as pre- 
dicted by Andersen et al. and that the results can be well 
described by the adaptation of the linearized Usadel equations 
to this new situation. However, there are important aspects of 
the results which are not explained on this basis. Firstly, there 
appears to be a longer period oscillation in T c , which is visi- 
ble in both data sets as an upturn in the trend in T c for N > 
5. Perhaps, more significantly, there appears to be a parity- 
dependence of the resistive transition width (AT C ) as shown 
in Fig. [5Jc), which suggests that the nature of the supercon- 
ducting transition is being affected. 





nm 


Vf I 
xlCfms- 1 K 


(T 

xlO 6 (fim) 


D 

- 1 xlO^nvV 1 


£ 

nm 


Nb 


5.8 


0.3 


1.9 


2.2 


2.2 


Fe 


3.7 


2.0 1000 


6.6 


18 


2.7 


Cr 


9.6 


2.0 


2.2 


6.0 


0.9 



TABLE I: Important parameters used/calculated in this letter. 
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